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The infrared features of Landau gauge QCD are studied by lattice simulation of /3 = 
6.0, 16 4 , 24 4 , 32 4 and (3 = 6.4, 32 4 , 48 4 . We adopt two definitions of the gauge field; (1) U 
linear and (2) log U, and measured the gluon propagator and ghost propagator. The infrared sin- 
gularity of the gluon propagator is less than that of the tree level result, but the gluon propagator 
at momentum remains finite. The infrared singularity of the ghost propagator is stronger than 
the tree level. The QCD running coupling measured by using the gluon propagator and the ghost 
propagator has a maximum a s (p) ~ 1 at around p — 0.5 GeV and decreases as p approaches 0. The 
data are analyzed with use of the formula of the principle of minimal sensitivity and the effective 
charge method, and by the contour-improved perturbation method, which suggest the necessity of 
the resummation of the perturbation series in the infrared region together with the existence of the 
infrared fixed point. The Kugo-Ojima parameter is about -0.8 in contrast with the theoretically 
expected value of -1. The color off-diagonal part of the ghost propagator in the Landau gauge is 
consistent with zero, and its fluctuation can be parametrized as a constant/(ga) 4 . 

PACS numbers: 12.38.Gc, 11.15,Ha, 11.15.Tk 



I. INTRODUCTION 

Two decades ago, Kugo and Ojima proposed a cri- 
terion for the absence of colored massless asymptoptic 
states in Landau gauge QCD using the Becchi-Rouet- 
Stora-Tyutin(BRST) symmetry They suggested to 
measure the two point function of the covariant deriva- 
tive of the ghost and the commutator of the antighost 
and gauge held 

Q 

= 1 Y^e-^-^itvtx^D^ -L- [A v ,X b ]) }.(1) 

x,y \ ' X V 

They showed that u ab {0) = S ab u(0) satisfies 



l + «(0) = 




where Zz is the gluon wave function renormalization fac- 
tor, Z\ is the gluon vertex renormalization factor, and Z3 
is the ghost wave function renormalization factor, respec- 
tively. Kugo claimed that confinement is realized cither 
by (1) Z x = and Z 3 =finite or (2) Z 1 /Z z = but(and) 
Z3 = 0. The divergence of Z3 implies u(0) = —1 and 
the presence of a long-range correlation between colored 
sources. 



As shown by Gribov 2], the Landau gauge is not 
unique and the uniqueness of the gauge field can be 
achieved via restriction to the fundamental modular re- 
gion (FMR), i.e. the region where the norm is the ab- 
solute minimum 0,13. We adopted the smearing gauge 
fixing to make the configuration close to the FMR. 
We observed proximity of the gauge configurations with 
and without smearing gauge fixing, but the overlap of the 
Gribov region and FMR remains to be investigated. 

The confinement scenario was recently reviewed in the 
framework of the renormalization group equation and 
dispersion relation j(| . It was shown that the gluon dress- 
ing function Za{q 2 ) defined from the gluon propagator of 
SU(n) 

D^{q) = £ e-^TriA^AM) 

X=x,t 

= (S^- q J^)D A (q 2 ), (3) 

as Z A {q 2 ) — q 2 D A (q 2 ), satisfies the superconvergence re- 
lation, and the gluon dressing function at zero momen- 
tum does not necessarily vanish as Gribov and Zwanziger 
conjectured, but it could be finite. A systematic study 
of lattice data indeed establishes the infrared finiteness 
of the gluon propagator Q • 

The ghost propagator is defined as the Fourier trans- 
form of an expectation value of the inverse Faddeev- 
Popov(FP) operator M. 

D a G b (x,y) = (triX^KMiU})- 1 ^}}, (4) 
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and whether it satisfies the superconvcrgcncc relation is 
not clear. In maximal abelian gauge, it is conjectured 
that the off-diagonal gluon and off-diagonal ghost sat- 
isfy the superconvergence relation, but the diagonal ghost 
does not, and it is the source of the long-range correla- 
tion. We remark that Nishijima proposed a sufficient 
condition for the confinement as Z% = 0, based on the 
convergence of the spectral function [8j . 

The non-perturbative color confinement mechanism 
was studied with the Dyson-Schwinger approach and 
lattice simulations E El 111 III 111 ES E3 • We pro- 
duced SU (3 ) ga uge configurations by using the heat-bath 
method 0,^3 and performed gauge fixing ^3- We an- 
alyzed lattice Landau gauge configurations of j3 = 6.0, 
16 4 ,24 4 ,32 4 and (3 = 6.4, 32 4 , 48 4 produced at KEK. 
Progress reports are presented in [l3j and an extensive 
report will be published elsewhere. The gauge field is 
defined from the link variables as 



log U type: U x ^ = e J 
U linear type: A XifJ , = 



/It — _ A 



trie 



The fundamental modular region of lattice size L is spec- 
ified by the global minimum along the gauge orbits, i.e., 

Al = {U\Fu(l) = Min g Fu{g)}, A L C tt L , where n L 
is called the Gribov region (local minima) and 
f2 l = {U\ - dD(U) > , dA(U) = 0}. 

Here Fjj(g) for the two options are fl9ll20| 



A9 t Ag 



(1 - ±Re tvUiJ 



• Mq) = Ex,, 

respectively. 

The Landau gauge fixing in the log U type is performed 
by Newton's method where the linear equation is solved 
up to third order of the gauge field, and then the Poisson 
equation is solved by the Green's function method for 
16 4 , 24 4 lattices and by the multigrid method for 32 4 , 48 4 
lattices EH- The gauge fixing in the U linear type is 
performed by the standard overrelaxation method. The 
accuracy of dA(U) = is 10~ 4 in the maximum norm. 

In the calculation of the ghost propagator, i.e. in- 
verse FP operator, we adopt the perturbative method 
with use of the multigrid Poisson solver |2l| , whose accu- 
racy was kept within 10~ 5 and we set 1% as an ending 
condition in the method [13. But later we also intro- 
duce the straightforward and preconditioned conjugate 
gradient (CG) methods [2^ for cross-checking of the cal- 
culation. In the preconditioned CG method, we take for 
the preconditioning operation the same truncated per- 
turbation series of inversion as that of the perturbative 
method. In the CG method, the accuracy of the conver- 
gence of the series is set to be less than 5% in the L2 
norm. 

We analyze these data using a method inspired by the 
principle of minimal sensitivity (PMS) and/or the effec- 
tive charge method l23l l24| and the contour- improved 
perturbation method |25j. 



In sec. II we explain the method of analysis, and in 
sec. Ill the lattice data are presented and compared with 
results of the theoretical analysis. We performed a cross- 
check of our program of SU(3) by performing the SU(2) 
lattice simulation. In sec. IV we show our data and com- 
pare with those results of other groups. Our conclusion 
and outlook are presented in sec. V. Some details of our 
method of calculating the FP inverse operator are given 
in the appendix. 



II. METHOD OF ANALYSIS 

In the infrared region, the QCD perturbation series 
does not converge and truncation of the renormalization 
group equation and resummation of the series to evaluate 
the renormalon effect was proposed 26]. On the other 
hand, the possibility of the presence of an infrared fixed 
point was discussed and methods to bridge infrared and 
ultraviolet regions via the renormalization group equa- 
tion were proposed [2! |24|. The method was recently 
applied to an analysis of lattice data [23 and succeeded 
in explaining qualitatively the data. We briefly review 
the method. 



A. PMS and the effective charge method 

In the PMS method, the nth-order approximation to 
the physical quantity 1Z is expressed by the correspond- 
ing series of coupling constant h^ n ' which is defined as a 
solution of 



0o 



x 2 p x (3 x 2 + fax 3 H h P n x n+2 



),(5) 



where the scheme-independent constant and logarithmic 
term are separated. 

When 1Z is the QCD running coupling from the triple 
gluon vertex from up to three-loop diagrams in the modi- 
fied minimal subtraction (MS) scheme, one sets the scale 
u 2 equal to the external scale q 2 and expresses 

Tl n = h {n) (I + Ax/i (n) + A 2 h (n)2 + ■■■ + A n h {n)n ) (6) 

516217 

where in the case of n = 3, A\ = 70/3, A2 = 

153 304676635 299961 81825 1— 1 

C3 ~ As ~ 6912 _ C3_ 64~ ~ C5 ~6i~ m 
When one defines yjj§{q) as a solution of 

V^/ms^) = ^0 log(g/A H g) 2 - -i log(j3 Vus(q)) ( 7 ) 

Po 

and expresses the solution of Eq. (J5J 

h(q) = 2/ms(<z){1 + J/ms^^/A) " OVA)) 2 ) 



+ vmii) 3 'o(pVpo - OVA)) 3 ) + ■••}, (8) 
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where (3 = 11, ^ = 102 , f3 2 



2857 



0s 



149753 



2 ' ^ 6 
3564^3, we can calculate 1Z via eq.JBJ. 

The parameter yj^(q) can be expressed as y defined 

as a solution of 



m a scries 



Tl(q 2 ) = B^q 2 ) + A n B n+1 (q 2 ) (15) 



R 1 " 1 



0o 



log(/3 y) 



and the function 



k{q 2 , y) = - + ^ log(/3 y) - 0o log(g 2 /A 2 
2/ Po 



MS' 



(9) 



(10) 



In |27| the parameter y is fixed via minimization of 
\(R n (y) - n 1 {y))/'R}{y)\ for each q 2 . There are subtle 
problems in fixing y of PMS in the low-energy region 
[2!j,We leave the fitting of the low energy region for the 
future and we fix y at ^ — 1.97 GeV by solving 



l/y = /3 log(M/A) 2 



Po 



(11) 



The choice of /i = 1.97 GeV corresponds to the inverse 
lattice unit 1/a of (3 — 6.0 and chosen by as the 
factorization scale of the effective charge method. When 
A = = 0.237 GeV we find the solution y = 0.01594, 
and we call this method of choosing y at a specific fi 

and define a s from ghost-ghost-gluon coupling, the MOM 
scheme. 



B. Contour-improved perturbation series 



Exact solution of the two-loop renormalization group 

q 2 

equation for x with variable t = log —tt 



0{x) 



dx 
~dt 



-x 2 {\ 



cx) 



^\og{q 2 /k 2 ) = --c\og{-+c). 

I XX 



(12) 



(13) 



The solution 
1 1 



x can be expressed as x(q 2 ) = 

, where W(z) is the Lambert W function 
c 1 + W{z) 

which satisfies W(z)e w ^ = z. We apply the dispersion 
relation and consider contributions on a cut of negative 
real axis in the space of q 2 , i.e. take q pure imaginary. In 
order to be consistent with the MS scheme, the variable 
z is defined as 



a (-l-M/2c) 



iKrr 



-V 

e v A 



-b/ c giKir 



MS 



where t 
-b/2c 2 



-Z(q 2 )e 



MS 



(2c/b)^ b A- 



MS' 



(14) 



K 



25| . The physical quantities 1Z are expressed 



2tt J_^c{l + W(Z(q 2 )e iKS ) 



III. LATTICE DATA 



A. Gluon propagator 



) n d6. (16) 



The gluon propagator on the lattice was measured by 
using cylindrical cut method |l2j . i.e., choosing momenta 
close to the diagonal direction. When the difference of 



their lattice constant 



1.885 GeV in = 6.0, 32 s 



64 and our a 1 — 1.97 GeV, 32 4 is taken into account, 
the data are consistent with |13(see Fig. P). 

The effective coupling y of the MOM scheme is calcu- 
lated from 



l/y = /3olog( M /A z ) 2 -^log(/?oy) 
Po 



(17) 



for [i = 1.97 GeV and A z = A^e 25085 / 37752 [2| ob- 
tained from the three-gluon vertex in Landau gauge per- 
turbation theory. The relevant solution of ea. (|17fl is 
y = 0.02227. 

The gluon dressing function is defined as Z A {q 2 ) = 
q 2 DA(q 2 )- Its inverse Z~ x is expressed in the two-loop 
perturbation series as [2£| 

A - V2 )- 13 / 22(1 „™,( 2 ) 

- ( 41245993 - ™ C 3)^ 2 ), (18) 
v 1874048 352 Sil " v ' 

where X z is a fitting parameter (see Fig. [2J. 

As shown in Fig. Q]and Fig. [31 the gluon propagators of 
24 4 , 32 4 and 48 4 as a function of the physical momentum 
agree quite well with one another and they can be fitted 
by 



Z-\q\y) 



D A {q 2 



Z(q 2 ,y)\ y =0 .o2227 _ Z A (q 2 ) 



r 



(19) 



in the q > 0.8 GeV region. At zero momentum, Da(0) 
decreases as the lattice size becomes larger. 

The gluon dressing function in the MOM scheme with 
y = 0.02227 fits the lattice data for q > 0.8 GeV, 
but there appears a discontinuity at z ~ 0.174 and 
z = 1/eMWe note that the q dependence of Z(q 2 ,y) 
in z < 0.17 is similar to 1/W (0.17 - z), that in 0.17 < 
z < 1/e is similar to Re[W_i(0.17 - z) + W-i(z - e -1 )], 
and that in z > 1/e is similar to 1/Wq(z — 1/e). 
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B. Ghost propagator 

The ghost dressing function is defined by the ghost 
propagator as G ab (q 2 ) = q 2 D G ab (q 2 ). In the MOM 
scheme, we fix the scale by choosing y as a solution of 

1/y = \og^/K gh f - ^ log(/3 y). (20) 

For n = 1.97 GcV and A gh = A^e 1757 / 2904 obtained 
from two-loop Landau gauge perturbation theory we find 
as a relevant solution y = 0.02142. 
The ghost dressing function is 



FIG. 1: The gluon propagator as the function of the mo- 
mentum g(GeV). /3 = 6.0, 24 4 (triangles) , 32 4 (diamonds) and 
13 = 6.4, 48 4 (stars) in logt/ definition. 



5271, 



,/ov,, 615512003 5697 >N 

v 7/incmo *7r</f ^• 3 ' 



7496192 



704 



1936 ' 



(21) 



100 
75 
50 
25 


-25 
-50 
-75 




0.1 



0.2 



0.3 



. 4 



0.5 



FIG. 2: The dressing function \ z Z(q 2 ,y) as a function of the 
variable z for a fixed y — 0.02227. The branches of the Lam- 
bert W funtion l/W (0.17-z), Re[W-i(0.17- z) + W-i(z - 



- 1 )], and 1/W (z 



are shown as dotted lines. 




FIG. 3: The gluon dressing function as the function of the 
momentum q(GeV). (i = 6.0, 24 4 (triangles), 32 (diamonds) 
and f3 = 6.4, 48 4 (stars) in the log U version. The fitted line 
is that of the MOM scheme. 



where A g is a fitting parameter. 

In Fig. H the 24 4 ,32 4 , and 48 4 lattice data are com- 
pared with 



Z g (q 2 , y)|j/=o.o2i42 



G(q 2 ) 



(22) 
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FIG. 



4: The ghost propagator as the function of the momen- 

4 oo4 , ] a ci A /lo4 



turn q(GeV). (5 = 6.0, 24 4 ,32 4 and (3 = 6.4, 48 4 in the logU 
version. The fitted line is that of the MOM scheme. Z s (q 2 ,y) 
is singular at A^jg- ~ 0.25 GeV which should be washed away 
by the nonperturbative effects. 



We observe that the agreement is good for q > 0.5 
GeV and better than the result of the PMS method of 
[23|. The ghost propagators were calculated by the per- 
turbative method and the straightforward and precondi- 
tioned CG methods. We found that the two CG methods 
are consistent and give better accuracy than the pertur- 
bative method in SU(2), and give correct result in the 
lowest momentum point of SU(3) 48 4 lattice. With the 
lowest momentum point of the 48 4 lattice calculated with 
the CG method, the whole data can be fitted by Eq. 
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C. QCD running coupling 

We measured the running coupling from the product of 
the gluon dressing function and the ghost dressing func- 
tion squared @- 

a s (q 2 ) = ^Z A {q 2 )G{q 2 f ~ (qa)-«">>+**°\ (23) 

The lattice size dependences of the exponents a d and 
ac are summarized in Table [I] 

TABLE I: The exponent of gluon dressing function near zero 
momentum an, near qa = 1 a' D , and the exponent of the 
ghost dressing function near zero momentum ag. log U type. 



p 


L 


an 


a! D 


OtG 


an + Zoic 


6.0 


32 


-0.375 


0.302 


0.174 


-0.03(10) 


6.4 


48 


-0.273 


0.288 


0.193 


0.11(10) 



The effective running coupling in the MS scheme is 
expressed by the series of coupling constants few as 

eq.© HQ. The result of the MOM scheme using 
y = 0.01594 is shown by the solid line in Fig. The 
lattice data of 24 4 , 32 4 and 48 4 and the MOM scheme 
agree in 0.5GeV< q <2GeV, but the fit is slightly over- 
estimated at q > 2 GeV. 
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FIG. 5: The running coupling a s (q) of (3 = 6.0, 24 4 (squares), 
32 4 (triangles), f3 — 6.4, 32 4 (diamonds) and 48 4 (stars) as a 
function of momentum (/(GeV) and the result of the PMS 

method in the MOM scheme. 

In the contour-improved perturbation series, the run- 
ning coupling in two loops is expressed as 

TZ 2 (q 2 ) = Bitf) + A 1 B 2 (q 2 ) + A 2 B 3 (q 2 ) + ■■■ . (24) 

The series truncated at order A\ is plotted in Fig. 
together with the 48 4 lattice data measured by using the 
log U definition. We observed that the fit is good for q > 
2GeV, but the non-perturbative effect is underestimated. 

Since in the perturbative calculation of the Landau 
gauge gluon vertex in the MS scheme, the A^jg is modi- 
fied to e 70 / 6/3 ° A^jg [33, we performed the same replace- 
ment. The result overestimates in q > 1 GeV region. 



Since the A3 term is not known and there are cancella- 
tions between sucessive terms, we fit the data by inclusion 
of half of the A 2 term. The result is shown by the dotted 
line in Fig. [fjj 
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FIG. 6: The running coupling a s (q) as a function of momen- 
tum q(GeV) of the /3 = 6.4, 48 4 lattice. The solid line is the 
result of 1Z 2 using A^jg- = 237 MeV. Dotted line is the result 
of e 70 ^ 6f3 °A-Ms an d including half of A%. The dashed line is 
the result of the MOM scheme. 

The similar nonperturbative effect was attributed to 
the gluon condensates in 0, |^| . The lattice data are 
qualitatively the same as the results of hypothetical r lep- 
ton decay [33, and the Dyson-Schwinger approach |33| . 

The lowest-momentum point of a s of (3 = 6.4, 48 4 be- 
comes consistent with results of other lattice sizes when 
it is calculated with the CG method. 



D. Kugo-Ojima parameter 

Our lattice data of (1) the Kugo-Ojima parameter 
c = — u(0), (2)the trace of the gauge field divided by 
the dimension e/d, and (3 ) th e deviation parameter h 
from the horizon condition [l3| are summarized in Table 

nn 

We observe that the Kugo-Ojima parameter of the U 
linear definition remains smaller than that of log U. The 
similar difference exists in the ghost propagator in the 
infrared region. 



TABLE II: The Kugo-Ojima parameter c in the U linear and 
log U versions. j3 = 6.0 and 6.4. 



p 


L 


ci ei/d 


hi 


C2 62/d 


h 2 


6.0 


16 


0.576(79) 0.860(1) 


-0.28 


0.628(94) 0.943(1) 


-0.32 


6.0 


24 


0.695(63) 0.861(1) 


-0.17 


0.774(76) 0.944(1) 


-0.17 


6.0 


32 


0.706(39) 0.862(1) 


-0.15 


0.777(46) 0.944(1) 


-0.16 


6.4 


32 


0.650(39) 0.883(1) 


-0.23 


0.700(42) 0.953(1) 


-0.25 


6.4 


48 






0.793(61) 0.954(1) 
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FIG. 7: The Kugo-Ojima parameter c as a function of 
logZ(1.97GeV). f3 = 6.4, 32 4 , 48 4 in log [/(stars), = 6.0, 
16 4 and 24 4 in log [/(triangles) and U linear (diamonds) ver- 



We plot in Fig. 0the value c as a function of log Z(fi 2 ) 
of (3 = 6.0, 16 4 , and 24 4 in log U and U linear definitions 
and (3 = 6.4, 32 4 and 48 4 in the log U definition. The 
value c of f3 = 6.0, 32 4 is almost the same as 24 4 . The 
value increases as the lattice size increases from 16 4 to 
24 4 and the extrapolation of the two definitions to those 
of a large lattice where c in log U and U linear seem to 
cross at c ~ 1. The linear extrapolation as the function 
of logZ(/i 2 ) is based on the factorizablity 



MOM' 



(25) 



when fi ~ 1.97 GeV, which allows us to express 



MOM' 



x (^"'(AmomV^^ms))' ( 26 ) 



) to its continuum 



cteSnitions will appear 



The difference of the speed of Zt(A- 
limit in the U linear and log U 
as a difference of the slope. However, the increase of c 
from 24 4 to 32 4 is small. The Kugo-Ojima parameter c 
of (3 = 6.4, 48 4 lattice calculated in the CG method is 
0.793(61), which is consistent to the result of the log?7 
definition of (3 = 6.0, 24 4 , and 32 4 lattice data. 



IV. SU(2) LATTICE DATA 

In the numerical simulation of the SU(2) Yang-Mills 
field, we took the U linear type gauge field and simulated 
(3 = 2.2 and 16 4 lattices. We took 67 samples taken after 
18 000 thermalization sweeps and up to 84 000 sweeps 
with intervals of 1000 sweeps |16|. To each sample, 
we performed parallel tempering gauge fixing (PT) and 
direct gauge fixing by the overrelaxation method (first 
copy). We define the scale by the relation 1/a = 0.938 
GeV and compare our data with those of [3 HH an d 



A. Gluon propagator 

The gluon propagator is shown in Fig.JSl We observe 
that above 1 GeV our data agree with |3J|, but in the 
infrared region our data have an enhancement. Suppres- 
sion at momentum is consistent with the data of l35l. 
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FIG. 8: The gluon propagator -Da(q) as a function of the 
momentum g(GeV) of PT samples. 



B. Ghost propagator 

The color diagonal component of the ghost propagator 
calculated in PT is about 6% less singular than that of 
first copy (Fig. 0). We performed the calculation of the 
FP inverse operator by using the CG method, since the 
matrix is symmetric positive definite. Our data in the in- 
frared are less singular than that of [23] • Although there 
are difference in the gauge fixing method (PT versus sim- 
ulated annealing), we do not understand the origin of the 
difference. 
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FIG. 9: The color diagonal ghost dressing function Do(q 2 ) X 
q 2 as a function of the momentum squared g 2 (GeV 2 ). The 
first copies (diamonds) are more singular than PT (triangles). 
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In the maximal abelian (MA) gauge, color symmetry is 
spontaneously broken by the ghost condensation [37ll38| . 
In the Landau gauge, there is no background field as in 
the MA gauge, but the structure of the color off-diagonal 
ghost propagator has not been known. In order to in- 
vestigate this problem, we measured color off-diagonal 
symmetric and antisymmetric (e a bcDQ(q,q)) matrix ele- 
ments, where D^(q, q) is the ghost propagator with color 
indices a and b. We observed that the color off-diagonal 
antisymmetric part is consistent with zero pointwise as 
is expected from the theoretical observation and that 
the color off-diagonal symmetric part multiplied by q A 
is consistent with zero over the ensemble average, but its 
standard deviation is almost constant in the whole mo- 
mentum region. The fluctuation can be parametrized as 
a/q 4 with a = 0.0176(28) GeV 2 , in the normalization tr 
t 2 = 1 . We observed the same qualitative features in the 
SU(3) 48 4 lattice, but a is about 1/9 of the SU(2) 16 4 
lattice, i.e., the fluctuation is statistical. 



C. QCD running coupling 

The result of the running is shown in Fig llOl As a re- 
sult of the difference in the ghost propagator, the running 
coupling is about 1/3 of |3J|- We observe suppression 
near momentum. 



2 

1 . 75 
1.5 
1 .25 
1 

. 75 
0.5 
0.25 



. 5 



1 . 5 



2.5 



FIG. 10: The running coupling a 3 (q) 
momentum g(GeV) of PT samples. 



as a function of the 



D. Kugo-Ojima parameter 

The Kugo-Ojima parameter c of the PT samples was 
0.690(52) and that of the first copy was 0.722(68). This 
difference is qualitatively the same as that of the ghost 
dressing function at momentum. 



V. CONCLUSION AND OUTLOOK 

There are two aspects of color confinement: i.e., (1) 
the presence of long-range correlation between colored 
sources and (2) the absence of massless gluon poles. The 
Kugo-Ojima criterion is a sufficient condition for the two 
aspects, but the lattice data do not verify that these cri- 
teria are satisfied. 

A new method of FMR gauge fixing in SU(2) is re- 
ported in ^(|. We observe that the gluon propagator is 
suppressed at zero momentum in SU(2), i.e., the expo- 
nent o>d < —0.5, in contrast to the SU(3) case where 
an > —0.5. In the simulation of SU(2), we observed dif- 
ferences in the Kugo-Ojima parameter of the configura- 
tion in the FMR and of copies randomly produced in the 
Gribov region. The Gribov copy affects the Kugo-Ojima 
parameter, and in the ghost propagator in the infrared 
region, the difference is about 4%. Color SU(3) contains 
I, U and V SU(2) spin components, and we expect that 
the Gribov ambiguity is the same order. 

In the lattice data, the singularity of the ghost propa- 
gator is stronger than the tree level and that of the gluon 
propagator is weaker than the tree level. The dependence 
on the U linear or log U definition of the gauge field is 
small in the gluon propagator consistent with |39j, but 
not negligible in the FP inverse operator. 

We aimed at detecting in the lattice dynamics a sig- 
nal of the Kugo-Ojima confinement criterion derived 
in the continuum theory, formulated with use of the 
FP Lagrangian and BRST symmetry. We also noted 
that Zwanziger's horizon condition, based on the lat- 
tice formulation, coincides with the Kugo-Ojima crite- 
rion pjl Il3|. However, our present data are not satisfac- 
tory to prove or disprove the confinement criterion. The 
color off-diagonal antisymmetric part of the ghost propa- 
gator HtJEI vanishes in the Landau gauge, but the off- 
diagonal symmetric part has fluctuations proportional to 
{qa)-\ 

Although there are problems in fixing y of the PMS 
in the low-energy region, an extension of the effective 
charge method is a possible solution. In an extension 
of the solution of the two-loop renormalization group 
equation expressed by the Lambert-VF function, a so- 
lution of Pade approximant of the three-loop renormal- 
ization group equation was shown (TO and numerical 
calculation was done for Nf > 3 42] . In the ana- 
lytical perturbation theory approach in one loop, one 
predicts |43| a non-perturbative infrared fixed point of 

a s (0) 1 1 ^ , 
= — — = ; . Extension to the two- loops 



47T 

is discussed in 



411 . For Nf — 0, one needs continuation. 



There is a conjecture, in combination with the confor- 
mal relation, that continuation from Nf in the confor- 
mal window(4 < Nf < 6) to Nf — would be possible 

We remark that the Orsay group analyzed QCD run- 
ning coupling in the Landau gauge from three gluon 
vertex. They separate the momentum space into q < 
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0.8GeV and 1.5GeV < q and fit the lower momentum 
region by the instanton liquid model and the higher mo- 
mentum region by including 1/q 2 power correction due 
to gluon condensates ^(|. We did not take a s (0) = 0. 
Agreement of the lattice results of QCD running coupling 
in q > 0.8 GeV region and the three-loop perturbation 
theory is reported in |47|. 

We observed that the contour-improved perturbation 
theory performs a resummation of the perturbation series 
and that we can understand qualitatively the Landau 
gauge lattice QCD data via these methods. 
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APPENDIX: THE NUMERICAL CALCULATION 
OF THE FADDEEV-POPOV INVERSE 

In this appendix we briefly explain the numerical 
method of calculating the Faddeev-Popov inverse. 

1. Perturbative method 

The ghost propagator, which is the Fourier transform 
of an expectation value of the inverse Faddeev-Popov op- 
erator M = -3D = -d 2 {I - M), 

D a G \x,y) = (tr(\'*x\(M[U])- 1 \\ b y)) ) (A.l) 

where the outmost (• ■ ■ } denotes average over sam- 
ples U, is evaluated as follows. We take the plane 
wave for the source bW = \ b e lqx and get the solu- 
tion = bW. Here 0W = (— A) _1 bW. We cal- 
culate iteratively <p [ ' l+1] = Mcj>M(x)(i = I,-- - ,k - 1). 
The iteration was continued until the maximum norm 
Max x |(/)[ fc ](x)|/Max x |X;tri 1( /' [il (x)| < 0.001 - 0.01. The 
number of iterations k is of the order of 60, in SU(2), 16 4 
lattice, and of the order of 100 in SU(3). We measure 
also the h norm |||<^(x)|||/||| YllZl 0^ (x)[[|. 

We define $ b (x) = J2i=i 4^ ( x ) an d evaluate 
{\ a e lqx , <J> b (x)) as the ghost propagator from color b to 
color a. 

In the low-momentum region of SU(2) we observed a 
specific color symmetry violation pattern, and in the case 



of SU(3) relatively large color off-diagonal matrix ele- 
ments suppressed the color diagonal matrix element. For 
a cross-check of the perturbative method, we adopted the 
straightforward conjugate gradient method and the pre- 
conditioned conjugate gradient method 22] in which the 
truncated perturbation series is used for the precondi- 
tioning. 



2. Preconditioned conjugate gradient method 

We define M. = —d 2 (I — M) and define the truncated 
M.^ 1 which is used in the perturbative method as B^ 1 = 

{I +M H h M m - 1 )(-A)- 1 . First we choose x° and 

define r° = b — A4x°. Using the multigrid Poisson solver 
we calculate the perturbation series 

r° = £? -1 r° (A.2) 

and define p° = f°. 

Then we begin the iteration for k = 0, 1, • • • , 

a k = -(r k ,r k )/(p k 7 Mp k ), (A.3) 



(A.4) 



r k+1 = v k + a k Mv k . (A.5) 

We check the norm of r fc+1 , and if it is not small, we 
calculate the perturbation series f fc+1 = B~ 1 r k+1 as be- 
fore. We define 

/3 fe = (f fc+1 ,r fe+1 )/(f fc ,r fc ), (A.6) 



p fc+i =f fc+i +/3fcpfc; (AJ) 

and go back to the beginning of the iteration cycle. By 
choosing a sufficiently large number of m, the conver- 
gence occurs after a few iteration cycles. 

The preconditioned method makes the li norm conver- 
gence faster than the straightforward conjugate gradient 
method, but its maximum norm is larger than that of 
straightforward method. The solution agrees with the 
straightforward conjugate gradient method within errors 
in the whole momentum region, but disagrees with the 
perturbative method in the lowest momentum point of 
(3 = 6.4, SU(3) 48 4 lattice. 
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